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We present results obtained in QCD with two flavors of non-perturbatively improved Wilson 
fermions at finite temperature on 16 3 x 8 and 24 3 x 10 lattices. We determine the transition 
temperature in the range of quark masses 0.6 < m x /m p < 0.8 at lattice spacing awO.l fm and 
extrapolate the transition temperature to the continuum and to the chiral limits. We also discuss 
the order of phase transition. 
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1. Introduction 

Knowledge of the QCD phase structure is important in order to understand the physics of 
the early Universe, structure of neutron stars, processes in heavy ion collisions, etc. In order to 
obtain predictions for the real world from lattice QCD, we have to extrapolate the lattice data to the 
continuum and to the chiral limit. Determination of the transition temperature T c was studied by 
many groups. T c from staggered fermions with 2 and 2+1 flavors of dynamical quarks [|l|, ||, ||, |J, 
and from Wilson fermions with 2 flavors [J|] was found to be consistent in the chiral limit. Edwards 
and Heller [[J determined T c for L t = 4 and 6 using nonperturbatively improved Wilson fermions. 
We compute T c on finer lattices with L t = 8 and 10 and extrapolate our high statistics data to the 
continuum limit. Our results for L t = 8 at /3 = 5.2 were previously reported in Ref. [J^j. 

2. Simulation 

We use the Wilson gauge action and the non-perturbatively improved Wilson fermions with 
c sw given in Ref. [^]. Configurations are generated on 16 3 x 8 (/3 = 5.2 and 5.25) and 24 3 x 10 
(j8 = 5.2) lattices at various K. The values of K and the corresponding number of configurations 
for 24 3 x 10 and 16 3 x 8 lattices can be found in Table 1. We use T=0 results obtained by QCDSF- 
UKQCD to fix the scale. The contour plot of lines of constant r^/a and m K /m p is shown in Fig. [I]. 
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Table 1: Simulation statistics. 



3. Critical temperature 

We use the Polyakov loop susceptibility to determine the transition point. We identify the 
critical value of K, K t , as a point, where the Polyakov loop susceptibility reaches its maximum. 
Applying a Gaussian fit in the vicinity of the maximum (as shown in Figs. ^, || and|4j) we find K t = 
0.13542(6) at /3 = 5.2 on 24 3 x 10, K t = 0.13444(6) at jS = 5.2 on 16 3 x 8, and K t = 0.13406(6) at 
j3 = 5.25 on 16 x 8. Using the Sommer scale ro=0.5 fm and interpolating r^ja to the critical point, 
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Figure 1: Lines of constant r^ja (solid lines) and constant m K /m p (dotted lines) at T = 0. Crosses and 
circles indicate points where simulations are done. 



we obtain for the critical temperatures T c = 197(2)MeV (at m K /m p = 0.63(2)), T c = 212(2)MeV 
(at m n /m p = 0.76(1)), and T c = 217(2)MeV (at m n /m p = 0.81(1)), respectively. 
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Figure 2: Polyakov loop (left) and its susceptibility (right) at /3 = 5.2 on 24 3 x 10 lattice. 

In Fig. |we compare our result for the critical temperature with those of Refs. [P, [||] and 
where we have assumed t/o = 425 MeV, and m v and m ps are vector meson mass and pseudoscalar 
mass obtained on the zero-temperature lattice. Our result is in quantitative agreement with the 
results of the Bielefeld group and CP-PACS. This is reassuring, as [|I|], [||] and [||] work at larger 
lattice spacing. 



4. Continuum and chiral limits 

In order to estimate T c in the continuum and the chiral limits we apply the extrapolation for- 
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Figure 3: Polyakov loop (left) and its susceptibility (right) at /3 = 5.2 on 16 3 x 8 lattice. 
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Figure 4: Polyakov loop (left) and its susceptibility (right) at /3 = 5.25 on 16 3 x 8 lattice. 



0.65 
0.6 

CJ 0.55 

o 

^ 0.5 
0.45 h 
0.4, 



▲ 


1 1 

this work (clover, L y =8) 


1 1 


♦ 


this work (clover, L T =10) 






E-H (clover, L T =6) 




□ 


Bielefeld (p4 fat improved, L T 


=4) 
















$ ^ * ^ 






I I 


1 1 



0.2 



0.4 



0.6 

(m /m ) 

v ps \' 



0.8 



0.25 
0.2 
,0.15 
0.1 
0.05 
0. 



i — ffl — 







^ this work (clover, L T =8) 
^ this work (clover, L T =10) 
E-H (clover, L T =6) 
□ Bielefeld (p4 fat improved, L T =4) 
CP-PACS(RG,clover, L T =4) 
J I l_ 



0.2 



0.4 



0.6 

(m /m ) 

v ps V 



O.X 



Figure 5: (left) The transition temperature T c in units of the string tension a(T = 0) as a function of m ps /ffi v 
and (right) T c in units of m v as a function of (m ps /m v ) 2 (right). The data is from this work (the up triangles 
and the diamond), Ref. [|l]] (the squares), Ref. [||] (the circle), and Ref. Q| (the down triangles). 
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mula, which takes into account the lowest order corrections coming from the finite lattice spacing 
and quark masses. We use two extrapolation functions: the first one is of the bare quark mass type 
(TYPE 1) and the second one is of the pseudoscalar mass type (TYPE 2). 

TYPE 1 : 7> = T c m ^°r + C a (-) 2 + C q (- - , (4. 1) 

TYPE 2 : T c r = (T c r ) m ^ + C a (-) 2 + C q (m ps )& . (4.2) 

where j 1 ™'""^' coiTesponds to the extrapolated value of the critical temperature, and j8 and 8 are 
the critical indices. Since 0(a) improved Wilson fermions are employed in this studies, we expect 
corrections of 0(a 2 ) only. 

We make an attempt to fit four values for T c tq (see Table 2), obtained at rather large quark 
masses, in order to estimate the extrapolation parameters used in Eqs. ( |4.1||i~2| ). If the transition in 
two-flavor QCD is of the second order, then the transition is expected to belong to the universality 
class of the 3D 0(4) spin model with 1//3 5«0.54. Current lattice calculations suggest that transi- 
tion in QCD for physical quark masses is not a true phase transition but a crossover [jlj, ^, ||, |], |5j, [!(]] . 
Since results obtained by CP-PACS indicated a good 0(4) scaling [j|], we extrapolate the value of 
the critical temperature using the value of 0.54 as 1//35. We rely on the assumption that the fi- 
nite temperature phase transition in two-flavor QCD is second order. The best fit parameters are 
presented in Table 3. 



T c r 


a/r 




mps 


L t 


J3 


0.501(5) 


0.200(2) 


0.040(4) 


1.58(6) 


10(DIK) 


5.2 


0.539(5) 


0.232(3) 


0.094(4) 


2.25(4) 


8(DIK) 


5.2 


0.551(5) 


0.227(2) 


0.123(4) 


2.61(4) 


8(DIK) 


5.25 


0.57(2) 


0.290(9) 


0.174(12) 


2.74(5) 


6(E-H) 


5.2 



Table 2: The available data for the critical temperature T c ro. 
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Table 3: The best fit parameters for the case if the transition is of second order. 

We get the critical temperature in the continuum and in the chiral limits, 
TYPE 1 fitting: T c = 173(3)MeV, TYPE2 fitting: T c = 166(3)MeV. 

The value of T c extrapolated by using the TYPE 1 fitting agrees with the values obtained by 
other groups, but the value from the TYPE 2 fitting is slightly smaller. Both fits should give the 
same result. The difference gives an indication of the systematic error. 

Finally, we mention the results of Ref. [Q] which support the first order of the transition. If the 
transition is indeed of the first order, then 1/(58=1. In order to test the sensitivity of T c to the order 
of the transition, we also performed a fit with 1//35=1. 
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Table 4: The best fit parameters for the case if the transition is of first order. 

TYPE 1 fitting: T c = 192(7)MeV, TYPE2 fitting: T c = 182(4)MeV. 

These values are rather large compared to the value obtained from Nf = 2 + 1 staggered 
fermions in Ref. The results of our fits do not discriminate between the first and the sec- 
ond order phase transitions because of rather large errors in the T c r$ values. We are continuing 
simulations on 24 3 x 10 lattice in order to get a better precision of the T c value for this lattice. 
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